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Abstract 

In this paper we study the phase of self-similar solutions to general Nonlinear 
Schrodinger equations. From this analysis we gain insight on the dynamics of non- 
trivial solutions and a deeper understanding of the way collective coordinate meth- 
ods work. We also find general evolution equations for the most relevant dynamical 
parameter wit) corresponding to the width of the solution. These equations are 
exact for self-similar solutions and provide a shortcut to find approximate evolution 
equations for the width of non-self-similar solutions similar to those of collective 
coordinate methods. 
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1 Introduction 



The nonlinear Schrodinger equation (NLSE) in its many versions is one of the 
most important models of mathematical physics, with applications to different 
fields such as plasma physics, nonlinear optics, water waves and biomolecular 
dynamics, to cite only a few cases. In many of those examples the equation 
appears as an asymptotic limit for the slowly varying envelope of a dispersive 
wave propagating in a nonlinear medium [1,2]. A new burst of interest on non- 
linear Schrodinger equations has been triggered by the experimental achieve- 
ment of Bose-Einstein condensation (BEC) using ultracold neutral bosonic 
gases [3,4]. 

Many different forms of the Nonlinear Schrodinger equation appear in practical 
applications but many of them are included in the following family of equations 



Preprint submitted to Elsevier Science 



8 February 2008 




Am + V(x, t)u + g(\u\ 2 ; t)u, 



(1) 



where u(x, t) is a complex function to be determined, x E f2 C W 1 , V(x, t) is a 
real function called the potential and g is a real function satisfying g(|/; t) — > 
0, y — > accounting for the nonlinearities. In many physical applications Q = 
]R n and w has at least finite L 2 and H 1 norms in the spatial variables because of 
its physical meaning of mass and energy of the solutions. Most of the results to 
be presented in this paper can be extended without problems to the so-called 
vector case, when u = (ui, ...,um), however, to keep the formalism as simple 
as possible we concentrate here on the scalar case. 

Many questions may be posed from the mathematical point of view on the 
properties of solutions of Eqs. (1). Results are available for local and global 
existence, asymptotic behavior, etc ... [2,5,6]. In some very restricted one di- 
mensional situations it is even possible to find the analytical expression of the 
solutions by the use of the so-called inverse scattering transform or equiva- 
lent methods [7]. Finally, the so-called moment method allows in some specific 
cases to obtain rigorous results for the evolution of integral quantities related 
to the solution of Eq. (1) [8,9,10]. 

However, very little is known (rigorously) on the dynamics of solutions which 
are asymptotically non-stationary (e.g. when either the conditions discussed in 
Refs. [5,6] do not hold or there are more complex non-autonomous situations, 
etc). 

For those situations it is customary in the applied sciences to use the so 
called time-dependent variational method, collective coordinates method or 
averaged-Lagrangian method [11,12] to obtain some insight on the behavior 
of the solutions. All of the names refer to the same idea, which consists on 
rewriting Eq. (1) as a variational problem on the basis of the Lagrangian 
density so that the problem of solving Eq. (1) is transformed into the one of 
finding u(x, t) such that the action 



be one extremum. Of course, this problem is in principle as complicated as 
that of solving the original nonlinear Schrodinger equation (1). The idea of the 
method of collective coordinates is to restrict the analysis to a specific family 
of functions, that is to find the extremum of (2) over a specific family of func- 
tions u(x, t) = <f(x, qi(t), qi((t)', t) with fixed f and letting qi(t), qi((t) 
be unknown functions containing the dynamics of ip. Once the test function 
is chosen one may explicitly compute an averaged Lagrangian 






(3) 
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and the problem becomes a finite-dimensional one for the unknowns q±, qx- 
Thus, the evolution equations for the parameters are obtained by solving the 
Euler-Lagrange equations 

d (dL\ dL 

The accuracy of the results depends crucially on the proper choice of the 
family of test functions ip. In practice, this choice requires some knowledge of 
the solution obtained either from numerical simulations, asymptotic behaviors, 
or on the basis of more or less heuristic reasoning. Also the type of solutions 
studied usually have a simple form to keep the complexity of the calculations 
under control. 



There is an immense body of references on the application of the collective 
coordinate method to many specific applications, normally from a heuristic 
perspective. Two panoramas for its application in different fields can be seen 
in [11] and [12]. 

One drawback of the method is that in the framework of Nonlinear Schrodinger 
systems it is usually possible to guess an approximate form for the amplitude 
of the solution, but the choice of the phase is more difficult to justify. Usually 
a simple form is chosen based on the experience of the linear Schrodinger 
equations as will be discussed later. 

In this paper we propose a way to choose the phase systematically by solving 
exactly the phase equation corresponding to self-similar solutions. We also 
discuss how to use this phase to construct simple evolution equations for the 
parameters which improve the variational equations. In this analysis it will be 
also clear why the usual form of the method of collective coordinates is able 
to provide reasonable results for the dynamics. 



2 Formulation of the problem 



We start our analysis from Eq. (1) and split u into amplitude and phase by 
writing u(x,t) = a(x,t)e t ^ x ' t \ After substitution in Eq. (1) we obtain the 
following coupled system of equations 



da 1 . , 

_ = _(Va) (V0)--aA0, 

9 ~dt = 2 Aa ~~ 2 ( V< ^ 2fl ~~ V ( x ' t " >a ~ 



(5a) 
(5b) 
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From now on we concentrate on non-propagating solutions of Eq. (1), i.e., those 
for which the "centrum" of the wavepacket defined as Xj(t) = J Xj(a(x, i)) 2 d n x 
satisfies Xj(t) = 0. Because of the Ehrenfest theorem [13,14], for any type 
of nonlinear term g and for cuadratic polynomial potentials, the evolution 
of this quantity can be decoupled exactly from the "internal" dynamics of 
the solution so that starting from solutions u(x,t) with Xj(t) ^ we may 
construct new solutions u(x — Xj(t),t) satisfying Xj(t) = as discussed in 
[14]. Thus for the above mentioned type of problems (which include many 
problems of practical interest) we may concentrate on the analysis of this type 
of solutions without loss of generality. For nonlinear Schrodinger equations 
with algebraic potentials or polinomial potentials of degree higher than two 
our analysis will be applicable to non-propagating solutions. 



2.1 The phase of self-similar radially symmetric solutions 

As discussed above, the key idea of the method of collective coordinates is to 
restrict the analysis to a particular set of functions ip(x,qi(t), ...,q M (t);t). In 
this paper we will call quasi-solutions to these functions which are not true so- 
lutions of the original equations (1) but in some sense provide some estimates 
for the dynamics provided the parameters qj(t) are chosen optimally and the 
choice of the set of trial functions is appropriate. In the case of Schrodinger 
equations, the time-dependent variational method proceeds usually by choos- 
ing quasi-solutions of the form 



Thus, when non-propagating symmetric systems are considered there appear 
usually three free parameters A(t), w(t), f3(i) to be determined (they have the 
physical meaning of amplitude, width and chirp of the solution) once the 
motion of the center has been decoupled. The consideration of non-symmetric 
systems over n spatial dimensions leads to 2n + 1 independent parameters and 
can be handled in a similar way. The specific choice for the phase /3(t)r 2 is 
believed to be a good one on the basis of what it happens in linear systems 
and previous experience with the equation. In fact, this term is chosen to 
be quadratic by analogy with the Optics of gaussian beams, ruled by linear 
Schrodinger equations. In that field it is well known that gaussian beams must 
include a cuadratic phase term to be exact solutions of the linear propagation 
equations. One of the results of this paper will be to provide an understanding 
of why this approximation works so well for the nonlinear case. 

The main point of our analysis is to study self-similar solutions of the form 



<p = A(t)p(r/w(t)y mr2 . 



(6) 



u(x,t) = A(t)p(r/w(t))e^ (r '* ) . 



(7) 
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being A(t), w(t), p and (j>(r,t) real functions. 



When self-similar solutions to Eq. (1) exist our Eq. (7) will provide exact 
solutions. When self-similar solutions do not exist we will understand Eq. 
(7) as an ansatz approximating the true solutions in a sense similar to that 
of the quasi-solutions of the method of collective coordinates. The point we 
will develop in what follows is that the phase can be found in closed form 
independently of the specific form of the amplitude p. 

One key feature of NLS equations given by Eq. (1) is the conservation of 
the L 2 -norm \\u\\2 = (J \u(x, t)\ 2 d n x) 1 ^ 2 under time evolution. For self-similar 
solutions and choosing the norm value to be equal to one (this can be done 
without loss of generality by scaling appropriately the nonlinear term) this 
conservation law becomes 

A 2 (t) [ (p(r/w)) 2 d n x = 1. (8) 

Let us choose p such that ||p|| 2 = 1 (this can be done without loss of generality 
and determines uniquely A(t)), then 

A 2 w n = 1, (9) 

thus 

A=~w. (10) 

The usual situation in applied sciences corresponds to functions for which not 
only the L 2 -norm is well defined (and conserved), but also the H 1 norm is well 
defined. In fact, the gradient term / \Wu\ 2 d n x is sometimes called the kinetic 
energy, which must be finite for physically meaningful solutions. Thus, we will 
look for self-similar solutions of the form (7) with finite H 1 and L 2 norms. 
Situations where blow-up occurs require an specific analysis near the blow-up 
point where / |Vw| 2 — > oo [2]. 

Substituting (7) into (5a) gives 



w 

Aip(r/w(t)) -A(t)Mr/w(t)) = 



A{t)p r {r/w)<t> r {r/w) - -A(t)p(r/w) 



i n — 1 
<p rr H <p r . 



, (11) 



We now define a new variable q = r/w(t) and then Eq. (11) becomes, after 
using (10) 

*- + + t)** = 2,6,0 S + T*)' (12) 
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It is remarkable that this is a linear second order equation whose general 
solution can be easily found. Defining v(q) = <j)/(2ww) we get 



v a + 



n-1 2p 



n p q 



+ — )v = - + ^q. 
p ) 2 p 



(13) 



The formal general solution of Eq. (13) is 

C2(t) 



v(q) 



+ 



for any arbitrary function C 2 (t). Thus, integrating we get 

ds 



(14) 



(j)(q, t) = 2ww 



Ci(t) + C 2 (*) 



?n— 1 



(15) 



This is the general form of the phase of any self-similar solution to Eq. (1) with 
amplitude profile given by p. However, not all of these solutions correspond 
to physically interesting ones. To understand why let us compute 



K(t) 



\X7u\ 2 d n x 



2 i 2 i2 
Pr + P <Pr 



d n x. 



(16) 



As discussed above this integral (which is related to the H 1 norm of u) must 
be finite. However, it is clear that 



K(t)> 

= S, 



P 2 4> 2 q d 



roc 

n q= p 

Jo 

c 2 (t) 



C2(t) 



c 2 (t) 



-.n-1 



p(qY 



f-'p{qf 2 



+ 



n ' l dq I dQ 
q\ dq + q n+1 p{qfdq 



(17) 



where S n denotes the surface of the n-dimensional sphere of unit radius and 
S n = jg„ dVt is its measure. The most singular term in Eq. 17 is the first one. 
Taking into account that J Rn p 2 d n q = 1 we can use the Schwartz inequality to 
prove the divergence of the first term of Eq. (17) 



X™ ^ ^) \Io q n 1 p(q) 2 ^ Is™ i 



q n - l p(qf 



>sl 



p2qn " 1 ^^ dq = C2{t)Sl 



dq -> oo. (18) 



Which proves the divergence of K(t) unless C2(t) = 0. The second term in 
(17) cannot compensate this divergence since this would require p 2 ~ l/q n 
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for q — > oo which would lead to non-normalized solutions. Using this fact and 
changing back to physical variable r we get finally for the phase 

0M) = £r 2 . (19) 

Eq. (19) is one of the main results of this paper: the phase of finite- energy 
self-similar solutions is cuadratic irrespective of the many possible variations 
of the problem ( nonlinear ity, potential, spatial dimensionality). This fact ex- 
plains why the choice of quadratic phases as in (6) in time- dependent vari- 
ational method leads to very good results. This fact also clarifies why these 
approximate methods lead usually to a prefactor (5 — w/ (2w), which is indeed 
the exact result. 

We must emphasize that phase choices such as (19) have been used previ- 
ously by other authors either as approximate expressions for the phase as 
comented before or as an specific choice leading to particular solutions (e.g. in 
the framework of the theory of self-focusing in NLS equations [2,16]). What 
we prove here is that all finite- energy self-similar solutions to Eq. (1) must 
have a cuadratic phase, which is a stronger result. 



2.2 Equations for the scaling parameter 



When self-similar solutions (or quasi-solutions in the sense discussed above) 
are considered it is possible to get closed equations for the evolution of the 
most relevant parameter w(t) which determines the width of the solutions. In 
this section we restrict ourselves to the autonomous case when V and g do 
not depend explicitly on t. Let us first consider the hamiltonian 

H=^J \Vu\ 2 d n x + J G(a 2 )d n x + J V(x)\u\ 2 d n x, (20) 

where G is a real function satisfying g = dG/d(a 2 ). It is easy to check that H 
is a conserved quantity which is assumed to be finite. For self-similar solutions 
of the form (7) and using again the definition q = r/w, we get 



H 



2 7r- 



dp 



dr 



dr 

,2 



d n x + f G(A 2 p 2 )d n x + A 2 f V{x)p 2 d n x 

JR n JR n 



+ / A 2 w n V (wq) p 2 d n q. 
Using the L 2 -norm constraint (9) we get 



(21) 
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H = 7 ± 1 f p 2 q d n q+-L [ p 2 4> 2 q d n q 
2w z JR™ H 2w z JR n H 

+w n [ G (p 2 /w n ) d n q+ [ V (wq) p 2 d n q. 

JR n ^ ' JR™ 

Finally, we use (15) and the fact that 



r roc r roo 

/ d n q = / q n ~ x dq / dtt = S n / q^dq, 

JK. n JO JS n JO 



[22) 



(23) 



to get 



H l r°° 



POO I POO 

/ P 2 q q n ~ l dq + -w 2 p 2 q n+l dq 
Jo 2 Jo 



S„ 2w 2 



+w n G (p(q) 2 /w n ) q n ~ x dq + £ V (wq) p 2 q n ~ x dq. (24) 
Differentiating (24) with respect to time we obtain 

n'ii POO POO POO . . 

- / PaQ^dq + wid / q n+1 p 2 dq + nw n ~ l w / G (p 2 /w n ) q n ~ l dq 
w 6 Jo Jo Jo 

poo nin n ~^i'n / \ f°° 

-w n ^^P 2 g(p 2 /w n )q n - 1 dq+ V'(wq)qwp 2 q n - 1 dq = 0, (25) 
JO w zn v ' Jo 



which can be written as a Newton-like equation 



f 

Jo 



q n+l p 2 J w 



(/o 00 Q n -y q dq) 



—nw 



w- 



n-l 



G [ — 

\w n . 



w n \w n . 



q n - x dq 



q n V'(wq)p 2 dq. 



(26a) 
(26b) 
(26c) 



The different terms in Eq. (26) have a very clear physical interpretation. The 
right-hand-side of Eq. (26a) is a repulsive term proportional to 1/w 3 which 
accounts for the tendency of the wave-packet to spread under the action of 
dispersion (this is mathematically described by the term with the Laplacian 
in Eq. (1). The second one (Eq. (26b)) is the self-interaction energy due to 
the nonlinearity. Finally the last term given by Eq. (26c) comes from the 
interaction with the external potential V. 

Eqs. (26) are of the same form as the equations arising in particular studies 
of Eq. (1) based on the method of collective coordinates of which Eqs. (26) 
are a generalization. From them it is possible to understand that when self- 
similar solutions exist what the variational method does is to fit somehow the 
right dynamics but using different values for the integrals depending on p. 
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Due to the lack of knowledge of the right self-similar profiles p, the integrals 
are estimated by using the ansatz tp. If the ansatz is constructed reasonably, 
e.g. taking into account the right asymptotic behavior of the solutions, we 
may expect only a quantitative discrepancy between the exact results for self- 
similar solutions and the estimates obtained from the collective coordinate 
method obtained from Eq. (26). 

The main advantage of our derivation of Eqs. (26) is that they are obtained 
systematically on the basis of the exact phase of self-similar solutions. These 
equations provide then a direct way to obtain estimates for the dynamics when 
self-similar solutions do not exist and allow to avoid the lengthy calculations 
of the method of collective coordinates. 

There is a lot of information contained in Eqs. (26). Specifically, they de- 
scribe many cases previously studied. For instance, the simplest case g(y) = 
goy 2 ,9o < 0, V = corresponds to the usual focusing cubic NLS. From there 
we obtain immediately the following condition for the collapse of self-similar 
solutions 

which allows to estimate the critical value of g for solutions with ||p||2 = 1. 
For instance taking gaussians one gets g c = —2n, which is the usual estimate 
obtained from time-dependent variational methods. Taking p to be the so- 
called Townes-soliton [2] leads to the exact result for the critical value of g 
[2]- 

Many other situations have been described in the literature. Taking V(r) = 
|0 2 r 2 and g(y) = goy 2 ,go G 1R we get the simplest situation of Bose-Einstein 
condensation in a symmetric trap. This problem was studied variationally in 
Ref. [17]. Our equations coincide with those for the spatially symmetric case 
discussed there. Also in Refs. [18,19,20] cubic-quintic nonlinearities where used 
to describe four-body collisions in Bose-Einstein condensates. Eqs. (26) give 
us the evolution equations found there by just choosing a nonlinearity of the 

type g(y) = g 3 y 2 + g$y A - 



3 Conclusions and discussion 



In this paper we have studied the phase of self-similar solutions to general 
Nonlinear Schrodinger equations. Our analysis provides a better understand- 
ing of the way collective coordinate methods work. Also we have used the 
results on the phase to find an effective equation describing exactly the dy- 
namics of self-similar solutions to the NLS equation and approximating the 
dynamics of other, more complicated solutions. 
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Our analysis on the phase can be easily extended to more complicated vec- 
tor systems. We expect that the procedure followed here may provide a basis 
for choosing systematically the phase for more complex ansatzs such as those 
which are not self-similar thus overcoming one of the difficult points of usual 
collective coordinate methods. It is worth mentioning that our analysis of the 
phase is applicable to non-autonomous situations such as the ones recently 
discussed in Refs. [21,22,23,24,25], which represent an emerging field of appli- 
cations of Nonlinear Schodinger Equations. 
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